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Abstract 


This paper examines the effects of magnetic field annihilation at 
a tangential or rotational discontinuity in a resistive plasma. The 
magnetic field intensity profile depends on : l) the field intensities , 

—Boo 3 fax from the current sheet (directional discontinuity) ; 2) the 
angle gj between and B^ ("Ti /'^ 00 =<1); and 3) the electrical 

resistivity, p . For a tangential discontinuity, the theory predicts 
a depression in B, centered at the discontinuity, when o> ><o 0 = cos - - 1 - 
3“B _J (2^1^ )-(2~ B x /^Bqjj) , and it predicts a monotonic transition 
fronTB^to when o> < o> 0 . The theory provides satisfactory fits 

to the magnetic field intensity and proton temperature profiles observed 
for two extremely broad D-sheets in the solar wind, using values of 
pt (t being the "diffusion" time) determined by the observed widths of 
the profiles. Assuming that t < 10 days, one obtains effective 
resistivities p>3 x 10 12 emu and p£ 2 x 1013 emu for these two D-sheets. 
Since most interplanetary "discontinuities" are at least an order of 
magnitude thinner than the D-sheet analyzed here, either the average 
resistivity at directional discontinuities is much lower than lO^ 2 emu 
or annihilation does not always occur at discontinuities. 



I. Introduction 


Burlaga and Ness ( 1968 ) and Burlaga (I 968 ) identified structures 
in the interplanetary magnetic field with dimensions .001 AU which are 
characterized by a region of depressed magnetic field intensity, B, 
centered at a discontinuous change in the direction of the field (see 
Figure 1). They presented evidence that these structures, called D-sheets, 
result from annihilation of components of the magnetic field at the 
directional discontinuity, and they suggested that the broad depression 
in B is due to dissipation of field components near this discontinuity 
resulting from a non-zero electrical resistivity of the plasma. This 
paper develops these ideas further, with emphasis on the effects of 
dissipation. 

Section II presents the theory for dissipation near a directional 
discontinuity and gives several theoretical profiles for B(x,t). Section 
III applies this theory to published observations of two D-sheets and 
obtains a lower limit on the effective resistivity at these D-sheets. Section IV 
shows that these resistivities are much larger than those calculated 
from the Spitzer formula, but comparable to the anomalous resistivities 
given by Scarf, 1970. 

The broad D-sheets discussed below are not frequently seen. 

Statistical studies have not been made, but it seems that they pass the 
earth at the rate of « l/month, approximately the same rate as shock 
waves. Thus the resistivities derived from these D-sheets are not 
representative of the interplanetary plasma; however, we believe that 
D-sheets, like shock waves, are of considerable physical interest in 
themselves, and it is in this spirit that this paper is written. 
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II. Theory 


The physical process which is formulated below is as follows. We 
assume that at a directional discontinuity in the solar wind, magnetic 
field is annihilated; the annihilation process determines the magnetic 
field intensity at the discontinuity. This intensity necessarily differs 
from the magnetic field intensity at jpo, and consequently there are 
gradients of B near the discontinuity. We postulate that the electrical 
resistivity near the discontinuity is non-zero and finite, but constant. 
Under these conditions, magnetic field is dissipated, and the D -sheet 
broadens with a characteristic profile. Vfe obtain a theoretical, profile 
which can be compared with observations and can be used to determine 
the effective resistivity. 

The geometry of the problem to be considered is illustrated in 
Figure 1. We consider a plane discontinuity surface corresponding to 
a directional discontinuity. Let the origin be on this surface, and let 

A * 

x be the direction normal to the surface. Let the magnetic field 
intensity and temperature depend only on x and t, i.e., B = B(x,t), 

T = ±T(x,t), where + refers to the region x>0 and - refers to the 
region x < 0, ^T is the sum of the proton temperature, ^T , and electron 

£T 

temperature, ^T e » Let B_ ss ®, where "*3^ are constant vectors along 

and let +B > B, thereby defining x. 

We consider D-sheets for which the total pressure is constant, 


l = ± nk ( ± T e + ) + ± B 2 /(8 tt ) = constant 


( 1 ) 



This is the case for tangential and rotational discontinuities. 

The momentum equation implies that the component of bulk velocity normal 
to the surface is constant. In a frame moving with the discontinuity 
surface, this is zero. The mass conservation equation then implies that 
dVdt = 0, where n is the number density. Thus n = n(x), independent of 
time. Changes in B as a function of time will be accompanied by changes 
in Tp + Te, in accordance with ( 1 ), assuming negligible heat conduction^ 
the magnetic energy which is dissipated appears as plasma thermal energy 
so that the total energy density does not change. Since the total 
pressure remains constant, there are no bulk motions. 

Assume that at the directional discontinuity a component of the 
magnetic field is annihilated as discussed by Burlaga and Ness ( 1968 ). 

For a tangential discontinuity, ~^B m and“E^ are parallel to the surface 
but not necessarily parallel to one another, and 



00 00 


where o> is the angle between +B and -B. For a rotational discontinuity, 
there is by definition a non-zero component of B normal to the surface, 
Ba, which is continuous across the surface ^assuming thermal isotropy 
(Burlaga, 1971a) J. Only the components of B parallel to the surface 
can be annihilated. Thus, for a rotational discontinuity, 
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( 3 ) 


B(O.t) “ ISaJ =| B* +( + Bn S m +V. + 2 1,1. • SiJ/li 


COS < 0 =~&, • + B, 

-lies — ||«o 


Both (2) and (3) can be expressed as B^ = ("b^ + + Bj/2; this relation 
does not involve the current sheet normal, which is not readily determinable 
for rotational discontinuities. The inner boundary condition is taken to 
be B (0,t) = Bjj^ "^ 00 /| + B j J for the solution + B(x,t) in x >0 and 
“B(°,t) = / | "B ffl | for the solution ~B(x,t) in x<0. Thus, we 

assume no change in magnitude across the current sheet, but we postulate 
a non-zero change in direction (see Figure 1). 

Assume that initially "^(x) and ~B(x) are constant, equal to 
2.co 3 respectively, and assume that the densities in the regions 

x > 0 and x < 0 are constant, equal to + n and "n, respectively. The 
densities do not change with time, as discussed above. 

The equation describing changes in the magnetic field is 


3£(x,t) a . 3B (x,t) 

at dx Sx ' 


(4) 


where 


"H = pAtt 


(5) 
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p being the electrical resistivity defined by J_ = E_/o in emu, Cowling 
(1957). This equation can be solved separately for the regions x > 0 
and x < 0, and since B ± is either constant or zero, (4) may be 
regarded as an equation for the intensity of the magnetic field 
component parallel to the surface, By(x,t). For a tangential 
discontinuity, this is simply the total intensity of the magnetic field 
| £|| | = ^(xjt). For a rotational discontinuity, the total, magnetic 
field intensity is *B(x,t) = (^B^ 2 + ^By 2 (xjt)) 1 / 2 , where is 

a constant not equal to zero, which is determined by the initial conditions. 

Solving (4) for constant T| subject to the above boundary conditions 
and initial conditions gives 


and 


+B II = B llmin + +A erf (C) x >° 


“ B it = By . + - A erf (r) x < 0 

'■ "rrn n N — 


where 


± k = 


- B. 


min 


( 6 ) 


(7) 


( 8 ) 


and 


C = x/(2VtvT) 


(9) 



Figure 2a shows B(£,t) for the case + B = ~B , B =0, for 

00 00 * Jl y 

several values of o> . Figure 2b shows the corresponding results for 
~B = + B /2, B x = 0. As illustrated by Figure 2a, there need not be a 
depression in B at x = 0 if "B^ ^ + B w . In fact, there will be no depres- 
sion whenever o> < (o = cos" 1 (3R/2 - (2R)-1), where R = "B / + B . For 
example <> if R = .5<>a > 0 - 105°o When = O^there is simply a gradual «, 
mono tonic increase from “B^ to +3^ * due to magnetic 

"diffusion". Figures 2a and 2b show the interesting facts that B(£) 
varies nearly linearly with £ for small £and 9 B/d£ ^0 at £ = 0. This 
can be shown by using the series expansion for erf (£), which gives 
for small £ 

^(C) = Bmin + ^ (2£MT (1 + £ 2 /3 +0 (£ 4 )) (10) 

In order to apply the above theory, it is necessary to determine 
whether the discontinuity under consideration is tangential or rotational. 
It is usually sufficient to determine whether or not the following 
necessary condition for rotational discontinuities is satisfied: 

\ - + 1 = 1 (-As. - A) (_i*£_) 1//2 ii - (ii) 

\ p. 4rr 1 j 

where “Py^ and Pj^ are the total thermal pressures parallel and 
perpendicular to “B, (Hudson, 1970) . If this condition is not satisfied, 
the discontinuity is tangential (assuming that discontinuities with 
shock-signatures have been excluded). Other bounds for rotational 
discontinuities under conditions at 1 AU are given by Hudson (1972). 
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III. Observations 


To relate the above theory to observations, one must consider 
that the discontinuity surface moves past the observer and that x can be 
in any direction. For a tangential discontinuity, the spatial scale of 
an observed profile is |x| = |v * xj r, where V is the solar wind velocity, 
t =|t - t Q |is the time relative to the time t G of the passage of the 
discontinuity, and x =-( B x B)/ | B x B |. 

Collecting the above results gives for a tangential discontinuity: 

±B(t) = Bmin + (=*£„ - B^jj) erf ( *3 t) ' (12) 

where 

% = * x/(2Vrjt') . (13) 

It is assumed that t > > <r» i.e. , that the elapsed time t since the 
formation of the current sheet, is larger than the time interval t during 
which the D-sheet moves past the spacecraft. If D-sheets, like directional 
discontinuities, are formed within .8 AU (Burlaga, 1971b) , then this 
condition is met by orders of magnitude. When (^ t)^ < < 3 one can 
use the series expansion for the error function (see Eq. 10) to relate 
% to the "half -width" =%• the ± t at which ±B = Brai n +( ± B - B^.J/2. 

The result is 

% = = M3/( Sr) a 1 *) 

4(±t) 

Comparing (13) and (l4) gives the important result that if T) 
is constant ( + T| = 1) ), and if t < < t, then + t^ = 
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Thus, ~B(t) and + B(t) are determined by four measurable constants: 

T , B min’ 3X1(1 ±B * In particular, if + B(t) is calculated from t w , 

B . . and ^B oo >then“B(T) is determined by only one additional number, 
min 00 

“B . Given ^B(t), the temperature can be calculated from (11) and 

GO 

compared with observations. 

April 3. 1966, Event . Figure 2 shows a D-sheet which was convected 
past Pioneer 6 on April 3> 1966. The data differ only slightly from 
the preliminary results in Burlaga (1968). As indicated by the bottom 
panel of Figure 2, the total pressure computed assuming a constant 
electron temperature, = 1.7 x 10^ °K, (Montgomery 1972; Scudder et 
al., 1972 ), is constant within the experimental uncertainties, as 
required by the theory in Section II. From Figure 2 we take B, 0 and 
cp equal to 4.4 7, -35 ° » l8o°, respectively for “B and 6.87, 20°, 350°, 
respectively for + B^ . In solar ecliptic (GSE) coordinates (x toward the 

A A 

sun, z normal to the ecliptic and y to form a right-handed system), 

B = (3.6, 0.0, 2.. 5) and ^ = (6.3*-l.l, + 2 , 3). The angle between 
"B and + B is o> = 175 0 + 15°. 

To determine whether or not the discontinuity is rotational., 
consider (11). If the pressure were isotropic, the change in the radial 
component of the solar wind velocity should be |"V - + V |= 115 km/sec 
if the discontinuity were rotational. In fact, the solar wind pressure 
is anisotropic, but the anisotropy was not measured for the event under 
consideration. Typically, however, 4rr (P|| - P x )/B 2 ^ .1 (e.g., see 
Hundhausen, 1972). Thus, if the discontinuity were rotational, one 
expects | V w - V *=* 105 km/sec. This is much larger than the observed 
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value , | AV | ^ 10 km/sec, so we conclude that the discontinuity is 
tangential. 

The theoretical value of B min , given by (2), is « 1.2/ for w = 175° 
B 00 ~ 4. 47, B» =6.87. The uncertainty in w is approximately 15° , 
and the corresponding uncertainty in is » +.37. Thus, the 

theoretical B min , computed from the annihilation hypotheses, is 
consistent with the observed B m j_ n = .9 7. 

^B(t) was computed from (12) using the observed B„_. _ and + B and 

nu- 11 00 

taking the 'half -width" of + B . to be t^ = 1.3 min. The result, shown 
in Figure 3, is in excellent agreement with the observations. Using 
only one additional quantity =4.47, ~B(t) was conputed from (12), 

with the result shown in Figure 3« Again, the result is in excellent 
agreement with observations. Note that no free parameters were used to 
obtain this result. 

The theoretical proton temperature profile (computed from the 
observed densities, the theoretical B(t), and the assumption that T g = 
constant = 1.7 x 10^ °K) is shown in Figure 3. There is qualitative 
agreement with observations in that the theory predicts a ma ximu m in 
proton temperature at the directional discontinuity, consistent with a 
large increase in temperature which is observed there. However, this 
conclusion is based on only one measurement, and the uncertainties in 
this measurement are very large. Furthermore, in the absence of electron 
measurements we cannot exclude the possibility that electrons are also 
heated at the current sheet. With such data, one cannot be certain that 
heat conduction is insignificant. 
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December 27, 1965 < event . Figure 4 shows a D-sheet which was 

convected past Pioneer 6 on December 27, 1965 . This was mentioned in 

the papers by Burlaga and Ness ( 1968 ) and Burlaga ( 1968 ), but the data 

here are more complete and are presented in an extended format of 

Figure 3. The total pressure (plasma plus field) is essentially constant 

if we assume T e = 1.50 x 10^ °K, + T e = 1.86 x 10^ °K, which are 
00 00 

both within the range of reported solar wind electron temperatures, 
(Montgomery 1972, Scudder et al. , 1972). For “B^and + B take, (B, 

0 se > cp sv P=(5.857, +10°, 100°) and (6.07, -20°, 185°), respectively. In 
GSE rectangular coordinates these yield "B = (-1.0, 5.67, 1.0)7 and 
+ B = (-5.61, -0.49, -2.0)7. For "V and + V we have taken (- 587 , -10, -45) 
km/sec and (- 563 , 20, -25) km/sec, respectively. To determine the nature 
of this discontinuity we use Eq. (ip.) neglecting the pressure anisotropy 
to compute the AV^. necessary for the discontinuity to be rotational: 
tSr = (51, 74, 35) km/sec. AV obseryed = (-24 + 5, -30 + 20, -20 + 20) 
km/sec. A pressure anisotropy of 1.4 x 10"^ ergs/cm-^ is required to 
force equality. This is clearly unreasonable since it exceeds the total 
plasma pressure for a (3 = 1 plasma with B = 5 . 87 , which is 1.35 x lO -1 ^ 
ergs/cnP. Even a mixture of more reasonable anisotropic effects plus 
experimental errors could not make the equivalence. We therefore 
tentatively conclude that this is a tangential discontinuity. The 
different electron temperatures adopted are also reasonably interpreted 
as being characteristic of two distinct regions of plasma with a non- 
propagating (tangential discontinuity) interface. As a check on this 
conclusion we test a necessary condition for a tangential discontinuity: 
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AV *(“B x + B ) = 0. For the values chosen — — —Ai. = -,09( 

- "• |avJ|Wb| 

which is consistent with our initial conclusion that the discontinuity 
is tangential. 

The components of B^- as interpolated across the discontinuity 
(Figure 4) give no new information about the classification, since by 
the generalization of Equations (2) and (3) the components of B^- are 
determined in the same way from 4 b for the two classes. The components 
of the vector B . are indicated by the d.ashed lines in the B field 
rectangular components panels of Figure 4. Since each point there 
represents < B^ >^_ over 15 secs, the agreement with theory is good. 

The total field plotted is < (£ B . ^ ) >. . Thus a reasonable empirical 

l ^ 

B min ( T = °) is (~ 2 ° 0 + 3*5» +0.5), yielding |B fflin |-4.1 + 0.27. Using 
this value of B^ yields a reasonable theoretical agreement for the 
remainder of the empirical data. 

We thus proceed on the supposition that the discontinuity is 
tangential with + A = 1.8/, "A = 1.657, , r ~ ~T = t = .5 rain, and 

a = .886, which yields the theoretical field profiles in the left 
hand panels of Figure 4. The theory qualitatively reproduces most 
temperature variations near and far from the sheet. 

The main discrepancies from theoretical profiles occur in the 
wings of the profile. A spatially dependent resistivity is probably 
required to explain the entire profile in detail. The basic features 
of this theory, suitably extended to accommodate p ( Q and/or number 
density) seems capable of explaining some of the observed D-sheet 
structures in the solar wind. 
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IV. Discussion 


It was shown above that for two D-sheets in the solar wind near 
1 AU, the observed B(t) profiles are quantitively in agreement with a 
theory which postulates magnetic field annihilation at a directional 
discontinuity and resistive dissipation of magnetic field gradients 
near the discontinuity. The theory predicts a depression in „B at the 
current sheet when ai > cos" 1 ( 3 R/2-(2R) -1 ) , where R = ~Bj + B. <1. 
Furthermore, it predicts that the "half-widths" of _ B(t) and + B(t) 
should be equal for constant |v • x| and T], ( + t^ = y = T ), and that 
t w is related to the electrical resistivity by (5), (13) and (14) i.e., 

p = 16 (V -x) 2 Tw 2 /t (15) 

The value of t, the length of time during which magnetic diffusion has 
been occurring, is unknown, but is probably less than the time required 
for an element of the solar wind to move from the sun to 1 AU which is 
presumably less than 10 days. If this is correct, the theory gives a 
lower limit on p 

p >16 (V * k) 2 Tw 2 /t 10 = (V • x) 2 t w 2 x 6.6 x 10 8 (16) 

where V is in km/ sec and t w is in min in the last equality. 

For the April 3, 1966, event: V = 4l5 + 5 km/sec, |v*x|= 128 km/sec 
which gives 


12 



p > 2 x 10^ emu 


(17) 


For the December 27, 1965, event, V = 575 + 15 km/sec, x = (0.32, 0.22, -.92), 
and t w = 0.5 min, giving 

12 

p > 3 x 10 emu ( 18 ) 

The values of p given above represent limits on the effective 
resisitivity derived from a macroscopic model with a very specific 
geometry. Currents do not appear explicity in the basic equations, 
but it can be shown that the constraint B = B(x)y implies that the 
induced electrical field, E, (E is not impressed externally) will not 
necessarily be along the current £, so the elementary concept of 
resistivity does not apply here. The situation is further complicated* 
in that u> e t e > > 1, where u) e is the electron gyrofrequency and t 
is the electron -proton collision time; Chapman and Cowling (1970) point 
out that "It is, questionable whether the concept of conductivity has 
any real significance when <u -t^is large". Finally, it is believed that 
in a collisionless plasma such as the solar wind the classical Spitzer- 
HSrm formula for a fully ionized plasma is probably not applicable because 
of the growth of various microinstabilities (Schindler, 1972 ; Scarf , 1970) . 
Comparing our empirical resistivities with those given by the Spitzer-Harm 
formula (Spitzer, 1962) we find that for T = 2 x lo5 % p*2 x 10 6 emu 
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and for T 6 x 10^ °K> p~ 10^ emu.; these values of p are 
on the order of six orders of magnitude smaller than the lower limi ts 
obtained from the observations. On the other hand, using the effective 
collision frequency for ion sound waves with strong turbulence (Scarf, 1970) 
we find if 2 x lO 3 ^ emu, Which is 1 to 2 orders of magnitude larger than, 
but still consistent with, our empirical lower bounds* 

We emphasize that the values of p given by (17) and (18) are not 
representative of the solar wind plasma near 1 AU. If they were 

representative, all directional discontinuities would show the effects 

described in Section II, contrary to observations. If annihilation and 

diffusion do occur at most directional discontinuities in the- solar wind, 

t w is typically < < 2 min. In fact, inspection of high resolution IMP I 

data shows that typically t w < 1 sec. The normal effective resistivity is 

correspondingly much smaller than the value 10 12 to lO 1 ^ emu derived 

above and could, in fact, be consistent with the Spitzer (1962)-Braginskii (1968) 

result. 

Finally, we note that the theory presented and applied in this 
paper may not be valid for all D-sheets. It was assumed that pressure 
is absolutely constant and that heat conduction is negligible, which 
implies no bulk motion and no change in density. However, observations 
of some D-sheets (Burlaga, 1970) show an appreciable enhancement in 
density, but no comparable increase in temperature at the corresponding 
discontinuities. This would occur if heat were conductedaway and material 
were conducted away and material were convected toward the discontinuity. 

The theory of this more complex process will be presented in a subsequent 
paper . 
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Figure 1 Shown is the geometry of the current annihilation model for 
either a rotational (ED) cr tangential (TD) discontinuity. 

If it represents a TD^B^x) = B(x) ; if a (ED) By (x) = B(x)-B^ 

The mag nitude of B(0) results from the annihilation of components 
of By after the suggestion of Burlaga and Ness ( 1968 ) . In 
general B(o) = B|| (0 ) + B x = l/2 ( + \) i this implies 
B|, ( 0) = 1/2 ( + B|| oo +"B ||w ) as is shown. 

Figure 2 Shown are theoretical profiles for annihilation 

models with no fluid motion relative to the surface and 
varying angles gj = cos”- L ( + ^ 11 * l^^for | B | B^ j = 1, 2 
on left and right hand sides respectively. £ is the dimensionless 
space-time variable defined in the text by Eq.. 9 . The profiles 
for rotational discontinuities may be determined via defining 

B x / B ll» and then |b(x) | = (B (j 2 (x) + B ^) 1 / 2 } giving a less steep 
profile. 

Figure 3 April 3, 1966, Event : Empirical field and particle profiles 
(dots) from Pioneer VI compared with theory (solid curves). 

Figure 4 December 27, 1965 , Event : Full data set used in determining 

type of discontinuity. Resultant comparison between relevant 
form of theory and empirical data are shown by solid curves. 
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FIGURE 1 
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